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Abstract 

We have performed the first-principles calculations to study the elasticity, electronic structure, 
and magnetism of InNCo3 and InNNi3. The independent elastic constants are derived from the 
second derivative of total energy as a function of strain, and the elastic modulus are predicted 
according to the Voigt-Reuss-Hill approximation. Our calculations show that the bulk modulus of 
InNCo3 is slightly larger than that of InNNi3 due to a smaller lattice constant for InNCo3. For 
InNCo3 the ferromagnetic state is energetically preferable to the paramagnetic state, while the 
ground state of InNNi3 is a stable paramagnetic (non-magnetic) state. This is due to the different 
strength of 2p-3d hybridization for the N-Co and N-Ni atoms in InNCo3 and InNNi3. 
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I. INTRODUCTION 



The ternary nitrides or carbides with the general formula AXM3 ( A: divalent or triva- 
lent element; X: carbon or nitrogen; and M: transition metal) are already known for several 
decades &H. These compounds crystalize in a cubic anti-perovskite structure (A: cube-corner 
position; X: body-center position; M: face-center position) and exhibit a wide range of inter- 
esting physical properties^, such as giant magneto-resistancei and nearly zero temperature 
coefficient resistivity!. They have renewedly attracted considerable attention due to the 
discovery of superconductivity at ~8 K in intermetallic compound MgCNi 3 i. 

Considering the Ni-rich composition, it is expected that the ferromagnetism could exist in 
MgCNi 3 . However, the absence of ferromagnetism was observed in experiment! for MgCNi 3 . 
From the electronic structures obtained by the first-principles calculationsB~0, the non- 
ferromagnetic ground state of MgCNi 3 is ascribed to a reduced Stoner factor that results from 
a strong hybridization between the Ni-3<i and C-2p electrons. For other Ni-based ternary 
carbides ACNi 3 (e.g., A = Al, Ca, In, Zn, and Cd), the first-principles calculations^^ also 
show that the ground states of these compounds are non-magnetic and the C-Ni bonding 
exhibits nearly same characteristics as the one in MgCNi 3 . Therefore, this indicates that the 
change of composition A could not induce the ferromagnetism in ACNi 3 . On the other side, 
it raises a question whether the change of composition X or M can lead to the appearance 
of ferromagnetism in AXM 3 or not. 

Very recently, the antiperovskite-type compounds InN^Cos and InN y Ni 3 (y ~ 1.0 and 
0.8, respectively) have been synthesized by solid-gas reactions of metal powders with NH 3 
and they have been reported to have spin-glass-like properties based on the measurements of 
temperature dependence magnetization.il! The recent first-principles calculations^! showed 
that the non-stoichiometry could affect the magnetic properties of ACNi 3 (e.g., AlCNi 3 and 
GaCNi 3 ) and suggested that the tendencies toward magnetism found in experiment sllHi! for 
these compounds should be explained by the deviation of the Ni/C atomic ratio from the 
ideal stoichiometry. To shed more light on the understanding on the magnetic properties 
of InNCo 3 and InNNi 3 reported in experimental, it is of great importance to theoretically 
study the electronic structures of these two compounds as well as the nature of the N-Co 
and N-Ni bondings. 

In order to completely understand the electronic structures and magnetic ground states of 
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I11NC03 and I11NM3 with cubic anti-perovskite structure, we carried out the first-principles 
calculations on these compounds using the pseudopotentials method with plane-wave basis 
set within the local density approximation and the generalized gradient approximation. Since 
the elastic properties of a solid are highly associated with various fundamental solid-state 
properties such as phonon spectra, specific heat, Debye temperature, and so on, we have 
calculated the independent elastic constant and the elastic moduli of InNCo3 and InNNi3. 



II. COMPUTATIONAL DETAILS 



All calculations on antiperovskite-type InNCo3 and InNNi3 were performed using the 
Quantum ESPRESSO codell, which is based on the density functional theory (DFT)0. 
The electronic exchange-correlation potential was calculated within the local density ap- 
proximation (LDA)Bi! and the generalized gradient approximation using the scheme of 
Perdew-Burke-Ernzerhof (PBE)@. The spin polarization was also considered in the calcula- 
tion in order to assess the magnetic properties of these compounds. Electron-ion interaction 
was represented by the norm-conserving optimizedil designed nonlocal pseudopotentials. 
The 4d electrons are explicitly included in the valence of In. The electronic wavefunctions 
were expanded by the plane waves up to a kinetic energy cutoff of 55 Ry. The k-point 
sampling in Brillouin zone (BZ) of simple cubic lattice was treated with the Monkhorst- 
Pack schemed and a 20x20x20 k-point mesh (i.e., 286 irreducible points in the first BZ) 
was used. The chosen plane-wave cutoff and number of k points were carefully checked to 
ensure that the total energy was converged to be better than 1 mRy/cell. The total energies 
are obtained as a function of volume and they are fitted with the Birch-Murnaghan 3rd- 
order equation of states (EoS)lzl to give the equilibrium lattice constant and other ground 
state properties. During the calculation of density of states (DOS), a dense k-point mesh 
of 30x30x30 is used, the total DOS is computed by the tetrahedron method^, and the 
atomic-projected DOS is calculated by the Lowdin populations^. 

For a cubic crystal, its independent elastic constants are Cn, Ci 2 , and c 44 . To determine 
the elastic constants of InNCo3 and InNNi 3 by means of the curvature of the internal energy 
versus the strain curves0'0, three strain modest are adopted and their nonzero strains 
are as follows: (1) e n = e 2 2 = <5, £33 = (1 + S)~ 2 ~ 1; (2) en = e 2 2 = e 33 = 5; and (3) 
e 12 = e21 = 5/2, e 33 = 5 2 /(4 - 5 2 ). The deformation magnitudes 5 from -0.012 to 0.012 
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in the step of 0.03 are applied in the first and second strain modes, and 8 from -0.04 to 
0.04 in the step 0.01 are adopted in the third strain mode. Once the independent elastic 
constants for single crystal properties are obtained through the above procedure, the elastic 
moduli (e.g., the shear modulus and the bulk modulus) of polycrystalline aggregates can be 
estimated according to the Voigt-Reuss-Hill approximation&Il. In the Voigt averageil, the 
shear modulus and the bulk modulus of cubic lattice are given by 

G v = \ [(en - c i2 ) + 3c 44 ] (1) 
5 

and 

B v = ~(cn + 2c 12 ), (2) 



while in the Reuss averageS they are given by 



Gr = -, (3) 

4(sn - S12) + 3s 44 

and 

Br 3s u + 6S12 ^ 
with the relations c 44 = s^ 1 , Cn — c 12 = (s u — s 12 ) _1 , and Cn + 2c 42 = (su + 2s 12 )~ 1 in the 

cubic lattice, where are the elastic compliance constants. Therefore, Gr and Br in the 

cubic lattice can be rewritten as 



Gr 

and 



4 j 3 

-(Cn - c 12 ) + -c 44 



(5) 



Br = - [( Cll + 2c 12 )] 

= By. (6) 

In the Hill empirical average@, the shear modulus and the bulk modulus are taken as 
G = \{Gy + Gr) and B = \{B V + B R ), respectively. Knowing G and B, the Young's 
modulus E and Poisson's ratio u, which are frequently measured for polycrystalline materials 
when investigating their hardness, can be calculated from the isotropic relations: 

and 

3B - 2G 

V ~ 2(3B + GY (8) 
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III. RESULTS AND DISCUSSIONS 



A. Structural properties 

In experiment with the powder X-ray diffraction patterns, Cao et allll have reported 
that InN y Co3 and InN^N^ (y ~ 1.0 and 0.8, respectively) have the cubic anti-perovskite 
structure with the space group 221(Pm3m) and the corresponding lattice parameters were 
3.854 A and 3.844 A, respectively. Starting from the experimental data, we have calculated 
the total energies of unit cell at a series of volumes for each compound in the paramagnetic 
(PM) and ferromagnetic (FM) states. The results are presented in Fig. |I| It is found that 
the energy difference between the FM and PM states is -0.0397 eV (-0.226 eV) in the LDA 
(GGA) calculations for InNCo3 and 0.0 eV for InNNi3. The total magnetic moment of 
InNCo3 is about 2.14 /i# (2.91 fis) and the local magnetic moment of each In ion is about 
0.69 fis (0.94 hb) in LDA (GGA) calculations. The total magnetic moment of InNNi 3 and 
the local magnetic moment of each Ni atom are zero. These indicate that the ferromagnetic 
state is energetically favorable to InNCo3 and the ground state of InNNi 3 is paramagnetic 
state (non- magnetic). The obtained equilibrium lattice constant (ao), bulk modulus (B), 
and first pressure derivative of bulk modulus (B') of InNCo3 and InNNi3 are listed in Table |. 
In our calculations, the predicted lattice constant of InNCo3 is slightly larger than that of 
InNNi 3 , which is oppsite to the trend reported in experimental. This may be due to the 
deviation of the Ni/N atomic ratio from the ideal stoichiometry in experiment for InN^Nis. 
In addition, it can be seen that the deviations of the LDA (GGA) lattice constants of both 
the InNCo3 and InNNi 3 with respect to the experimental values are less than 2.6% (1.0%). 
That is to say, the calculated equilibrium lattice constants of InNCo3 and InNNi3 are in 
excellent agreement with the experimental dataEl. 

B. Elastic properties 

The calculated independent elastic constants for single crystal of InNCo3 and InNNi3 
are listed in Table [TT]. Based on the Voigt-Reuss-Hill approximation&ll, the elastic moduli 
of InNCo3 and InNNi 3 are estimated and the results are listed in Table [TJ. For the bulk 
moduli of InNCo3 and InNNi 3 , the estimations based on the independent elastic constants 
agree well those obtained by the fit of the Birch-Murnaghan 3rd-order EoS. To the best of 
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our knowledge, no experimental data or theoretical results for the elasticity of InNCo3 and 
InNNi3 compounds have been reported up to now. Considering that the elastic properties of 
ZnNNi 3 , InNSc3, and InCNi 3 with cubic anti-perovskite structure have been studied recently 
and the theoretical results are available in literatureBBS, it will be meaningful to compare 
them with those of InNCo3 and InNNi 3 compounds. The order of bulk moduli of these 
five compounds from low to high is: 5(InNSc 3 ) < 5(InNNi 3 ) < 5(InCNi 3 ) < 5(InNCo 3 ) 

< 5(ZnNNi3). This could be understood from the trend in lattice constants (a) of these 
compounds (i.e., a(InNSc 3 ) > a(InNNi 3 ) > a(InCNi 3 ) > a(InNCo 3 ) > a(ZnNNi 3 )) as well 
as the relationship between bulk modulus and equilibrium volume (i.e., B ~ V The 
GGA lattice constants of InNSc 3 , ZnNNi 3 , and InCNi 3 are 4.411 A§, 3.77 A0, and 3.880 
A0, respectively. The order of shear modulus from low to high is: G(ZnNNi 3 and InNNi 3 ) 

< G(InCNi 3 ) < G(InNCo 3 ). Pugh0 has proposed that a high B/G ratio may be associated 
with better ductility whereas a low value would correspond to a more brittleness, and the 
critical value separating ductile and brittle materials is around 1.75. From the results of 
B/G ratio for InNCo 3 , InNNi 3 , InNSc 3 , ZnNNi 3 , and InCNi 3 , it is found that only InNSc 3 
can be classified as brittle materials and others may be ductile materials. Furthermore, 
InNNi 3 seems to be more ductile than InNCo3. 

For a cubic crystal, its mechanical stability requires that its three independent elastic 
constants should satisfy the following relations^: 

(cu - cia) > 0, en > 0, c 44 > 0, (en + 2ci 2 ) > 0. (9) 

These conditions also lead to a restriction on the magnitude of B: 

c l2 <B< c n . (10) 

The predicted values (see Table ||) for InNCo3 and InNNi 3 satisfy these conditions, 
indicating that cubic antiperovskite-type compounds InNCo 3 and InNNi 3 are mechanically 
stable. 

C. Electronic structures 

In order to understand the different magnetic ground states for InNCo3 and InNNi 3 , we 
examined the electronic structures of these two compounds. For the simplicity in discussion, 
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only the GGA results are presented below. The calculated electronic band structures along 
the high symmetry directions in the Brillouin zone are shown in Fig. 0. As discussed above, 
for InNCo3 the ferromagnetic state is energetically preferable to the paramagnetic state, 
and hence it is clearly seen that the spin-splitting occurs in the bands around the Ep. For 
InNCo3 the profile of majority spin bands looks roughly similar to the one of minority spin 
bands, however only one band in the majority spin bands crosses the Fermi level and the 
corresponding band with same dispersion in the minority spin bands is unoccupied. In addi- 
tion, two bands in the minority spin bands across the Fermi level of InNCo3 (see Fig. 0(a)). 
For InNNi3 the ferromagnetic state is not energetically preferable to the paramagnetic state, 
consequently, the majority and minority spin bands are degenerated, that is to say, no spin- 
splitting occurs in the band structure (see Fig. @(b)). It is interesting to note that the whole 
feature of majority spin bands of InNNi3 is very similar to the one of InNCo3. In order 
to reveal the detailed character of band structure, the total density of states (DOS) and 
the angular-momentum-projected DOS of each atom in InNCo3 and InNNi3 are presented 
in Fig. 0. For both InNCo3 and InNNi3, the four bands from -9 eV to -5 eV come mainly 
from the N 2p states and In 5s states, and the five bands roughly from -5 eV to -2.5 eV 
have significant contribution from 3d-t2 9 states of transition metal atoms (Co/Ni) and the 
5p state of In atom. For the bands from -2.5 eV to eV (i.e., the Fermi level), they are 
dominated by the 3d states of transition metal atoms (Co/Ni) and have small contribution 
from the 2p states of N atom. Because the number of 3d electrons in Co is one less than 
that of Ni, two minority spin bands composed of the 3d-t2 9 states around the Fermi level 
are unoccupied in InNCo3, while the counterpart in the InNNi3 are occupied. This also 
results in different behavior of spin-splitting in InNCo3 and InNNi 3 . Due to the significant 
spin-plitting around the Fermi level in InNCo3, the hybridization between the Co 3d and 
N 2p states in InNCo3 are slightly weaker than the one between Ni 3 c? and N 2p states in 
InNNi 3 (see Fig. |). 

The contributions of each kind of atoms to the DOS at the Ep of InNCo 3 and InNNi3 



are listed in Table [II. The total DOS at the Ep of InNCo3 is about 2.761 states/eV per 
formula unit (f.u.) in the GGA calculations, which is larger than that of InNNi 3 , and its 
main contribution comes from Co 3d states which accounts for 87%. For InNNi 3 in the GGA 
calculations, the contribution of Ni 3d states to the the total DOS at the Ep (i.e., 1.803 
states/eV.f.u.) accounts for 72%. These indicate the 3d states of transition metal atoms in 
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I11NC03 and I11NM3 play dominant roles in the total dnesity of states of these compounds. 

In order to understand the bonding nature among the ions in InNCo3 and InNNi3, we 
analyzed the charge density contours of InNCo3 and InNNi 3 in the(llO) plane, as shown 
in Fig. U From Fig. |], it is found that a certain amount of charges are accumulated in 
the intermediate region between N and Co atoms in InNCo3, and slightly more charges are 
accumulated intermediate region between N and Ni atoms in InNNi3 . This gives an evidence 
for the strong hybridization between N and transition metal (Co/Ni) atoms, indicating that 
the N-Co and N-Ni bondings exhibit strong covalent characteristics and the latter is slightly 
stronger than the former. The similar bonding characteristics for Ni-N atoms or Ni-C atoms 
were also found in other Ni-based ternary nitrides or carbides AXNi 3 0'0ilii. Therefore, 
our results suggest that the magnetic properties of InNNi 3 reported in experimentEl are very 
likely due to the non-stoichiometry effect, which was also found in the cases of AICM3 and 



IV. CONCLUSIONS 



In summary, we performed the first-principles calculations to study the elastic and elec- 
tronic properties of cubic antiperovskites InNCo3 and InNNi 3 . Based on the Voigt, Reuss 
and Hill bounds, the shear, Young's moduli and Poisson's ratio have also been estimated 
for the InNCo3 and InNNi3 polycrystals. The theoretically predicted equilibrium lattice 
parameters are in good agreement with the available experimental data. Our calculations 
show that the 3d states of transition metal atoms in InNCo3 and InNNi 3 play dominant roles 
near the Fermi levels. InNCo3 energetically prefers to the ferromagnetic state. The mag- 
netic ground state of InNNi 3 , which is same to other Ni-based ternary nitrides or carbides 
with a cubic anti-perovskite structure, is a stable paramagnetic (non-magnetic) state. This 
could be understood from that the hybridization between Ni-3<i and N-2p states in InNNi3 
is slightly stronger than the one between Co-3<i and N-2p states in InNCo3 because of the 
more 3d electrons in Ni. 
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TABLE I: Calculated lattice constants (a, in A), bulk modulus (B, in GPa), and the first pressure 
derivative B' of bulk modulus for InNCo 3 and InNNi 3 . AE™~ PM (in eV per formula unit)is the 
difference between the total energies of ferromagnetic (FM) and paramagnetic (PM) states for 
I11NC03 and InNNi3. The available experimental values are also listed. 

InNCo 3 InNNi 3 



LDA GGA Expt.0 LDA GGA Expt@ 

PM FM PM FM PM/FM PM/FM 

a 3.744 3.753 3.835 3.855 3.8541 3.784 3.882 3.8445 

B 255.78 243.06 211.29 194.17 226.91 179.93 

B' 4.497 4.483 5.570 5.568 4.761 4.281 

AE™~ PM - -0.0397 - -0.226 0.0 0.0 
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TABLE II: Calculated elastic constants (en, C12, and C44, in GPa), shear modulus (G, in GPa), 
Young's modulus (E, in GPa), and Poisson's ratio (y) of InNCo3 and InNNi3. The Voigt shear 
modulus {Gy, in GPa) and the Reuss shear modulus (Gr, in GPa) are also presented. For compar- 
ison, the elastic properties of InNSc3, ZnNNi3, and InCNi3 compounds with a cubic anti-perovskite 
structure are listed. 



Compound Method 


Cll 


C12 


C44 


Gy 


Gr 


B 


G 


E 


V 


B/G 


InNCo 3 


LDA 


389.11 


171.12 


102.55 


105.13 


105.03 


243.78 


105.08 


275.63 


0.311 


2.320 




GGA 


317.54 


126.76 


94.98 


95.14 


95.14 


190.35 


95.14 


224.67 


0.286 


2.001 


InNNi 3 


LDA 


356.77 


164.23 


69.06 


80.35 


78.12 


228.41 


79.24 


212.94 


0.344 


2.895 




GGA 


274.08 


131.20 


60.01 


64.58 


64.11 


178.83 


64.35 


172.37 


0.339 


2.779 


InNSc 3 


GGA§ 


238.57 


54.28 


90.76 


91.31 


91.31 


115.71 


91.31 


216.88 


0.188 


1.267 


ZnNNi 3 


GGA0 


354.28 


134.01 


48.06 


72.89 


62.05 


207.43 


67.47 


182.61 


0.353 


3.074 




GGA0 


364.20 


124.90 


32.69 


67.47 


46.09 


204.67 


56.78 


155.92 


0.373 


3.604 


InCNi 3 


LDA0 


414.67 


135.72 


68.02 


96.60 


85.56 


228.71 


91.08 


241.22 


0.324 


2.511 




GGA0 


344.20 


106.30 


62.68 


85.19 


77.31 


185.60 


81.25 


212.70 


0.309 


2.284 
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TABLE III: Total and partial density of states at the Fermi level (N(Ep), in states/eV.f.u.) for 
InNCo 3 and InNNi 3 . 

Total density Partial density of states 
of states In N Co 3 /Ni 3 
InNCo 3 LDA 3.811 0.073 0.337 3.401 
GGA 2.761 0.144 0.169 2.448 
InNNi 3 LDA 1.580 0.160 0.212 1.208 
GGA 1.803 0.180 0.256 1.367 
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Figure Captions 



• FIG. P Total energy (in eV per formula unit) versus the atomic volume (in A 3 ) for 
InNCo3 (upper panel) InNNi 3 (lower panel). 

• FIG. 0: (Color online) Electronic band structures obtained with GGA for the majority 
(red solid line) and minority (blue dashed line) spins of (a) InNCo3 and (b) InNNi 3 . 

• FIG. |3]: Total and partial density of states (DOS) of InNCo3 (left panel) and InNNi3 
(right panel) obtained with the spin-polarized GGA calculations. 

• FIG. [|: Charge density contours of the (110) plane for (a) InNCo3 and (b) InNNi3 
obtained with the spin-polarized GGA calculations. 
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